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A phenomenological theory of the incommensurate phase transition in thin films is presented for 
the case of the single-component order parameter. The order parameter distribution over the film 
is obtained under the assumption of slowly-varying amplitude. The size effect on the temperature 
of the phase transition into the incommensurate phase has been considered. 

I. INTRODUCTION 

Numerous recent experimental studies of thin films and small particles manufactured on the base of ferroelectric 
materials, such as BaTiC>3, PbTiC>3 and some others (see, for example, [1-8] and references therein) have shown that 
temperatures of phase transitions into a polarized state change noticeably (mostly, being reduced) with particle size 
or film thickness decrease, the phase transitions smeare out and the ferroelectric state is normally destroyed at some 
particular linear size. Spontaneous polarization, susceptibility, coercive and switching fields are also affected. Most 
of these features could be explained within theoretical models developed recently for ferroelectrics (see [9-15]). In 
particular, a phenomenological approach to the ferroelectric phase transitions in thin films and small particles was 
described in [9-11], while the Ising model was used in [12, 13]. The antiferroelectric phase transition in small particles 
and thin films was considered later in [17] on the base of the Landau theory. In all papers cited above, the polar states 
in relevant bulk materials were nonmodulated and the spatial variation of the spontaneous polarization and other 
parameters was induced by surface only (disregarding fluctuations). It is of interest to extend the phenomenological 
description of phase transitions in low-dimensional systems to transitions into modulated states; in particular, into the 
incommensurate phase. Note that some experimental studies of incommensurate thin films were made recently in [18]. 
The main difficulty of the phenomenological consideration of the size effect on the incommensurate phase transition 
is due to the fact that the incommensurate phase in the bulk is described by a nonlinear fourth-order differential 
equation for the order parameter contrary to the case of ferroelectrics or antiferroelectrics [19]. Thus, the boundary 
conditions for the order parameter, which were first suggested by de Gennes [20] and then were generally used in 
phenomenological approaches to the size effect, do not suffice. However, as it will be shown below, some assumptions 
allow us to reduce the problem to a simpler differential equation and to analyze the incommensurate phase transition 
quantitatively. The boundary conditions for this approximate equation are deduced from the conditions for the 
' fourth-order one. 

, In the present paper we suggest for the first time a phenomenological model of the transition into the incommensurate 
l 1 "^ ' phase in thin films. 
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II. THEORETICAL MODEL 
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Let us consider a thin film of thickness L. We assume that the x axis of the coordinate system is perpendicular 
to the film surfaces and the origin x = coincides with the film middle. Notice that due to the geometry chosen, 
the order parameter depends only on x. The order parameter r] can be the spontaneous polarization if the sample 
undergoes a transition into the ferroelectric phase at lower temperatures, but the further treatment does not depend 
significantly on its physical nature and may be applied to incommensurate phases of any physical nature. The 
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functional determining the equilibrium equation and the boundary condition can be written in the form 
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where Aq, B,C, D are phenomenological parameters, A$ > 0, B > 0, C < 0, D > 0, To is the characteristic tem- 
perature; the functions ijj(r](L/2)) and <fi(r)(—L/2)) determine the boundary conditions [22]. Having in mind that 
the symmetric solutions give usually the lowest magnitude of the thermodynamic potential in the case of symmetric 
samples, we consider below only even functions of x for the order parameter distribution: r}{—x) — tj(x) and only the 
region of x > will be considered explicitly. Varying the functional 0J, we obtain the Euler equation [21]: 
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describing the inhomogencous equilibrium state and the natural boundary conditions [21] 
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where F(r)(x),r)' ,r)") is the integrand of the functional JIJ, rf = -r^, rj' = We choose the simplest relevant form 
for ip: 
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Using the boundary conditions© and JSJ, and the functional Q), we obtain: 
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These boundary conditions can be used for numerical solving of the fourth-order equation J5J). However, we consider 
an approximate reduction of this equation to the second-order one, which will be solved exactly. 



III. SLOWLY- VARYING AMPLITUDE 



It is known that the order parameter in the bulk incommensurate phase oscillates with a steady amplitude according 
to the equation [19] r\ = ucos(qx), where the oscillation wave vector q is given by the relation 
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When the oscillation period in the bulk is implied to be much shorter than the film thickness, i. e. \D/C\ <C L 2 , we 
can assume that the restricted geometry does not influences noticeably the period of the order parameter modulation 
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and yields only alterations in the amplitude u(x) as a function of x. Then the x dependence of r\ in the film can be 
written as 



rj(x) = u(x)cos(qx). (9) 

Under the assumption of \D/C\ <C L 2 it is reasonable to assume in the spirit of the slowly- varying amplitude approach 
that variation of u(x) is smooth enough and its derivatives of order higher than the second one can be omitted. Then 
the relationship (JSJ leads to the following second-order differential equation for the amplitude u: 

d 2 u SB , C 2 - AAD , . 

dx^ = -SC U + ^CD~ U - (10) 

This equation along with relevant boundary conditions yields a boundary value problem for the amplitude u within 
the thin film. The boundary conditions for u can be deduced from the conditions for the fourth-order equation. 



IV. BOUNDARY CONDITIONS FOR THE SECOND-ORDER EQUATION 

The equation (|10|l must be supplemented with boundary conditions that can be derived from © and Q using 
the slowly-varying amplitude approach. We substitute the expression into these conditions and neglecting higher 
derivatives of u(x) we obtain a set of equations: 

<-l =q{C + Dq 2 )suJ q 4), (11) 



(C + 3Dq z ) — + a 
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^cosf^=-2^W<f). (12) 
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Eliminating the sine and cosine functions from (jl 1|> and ljl2|l and introducing the following notation for the logarithmic 
derivative, 

= 7 ' (13) 

we obtain a quadratic equation for the logarithmic derivative: 

a 2 

(C + Dq 2 ) ^-S 2 - a8 + (C + 3Dq 2 ) = 0. (14) 

Its solution reads 



y _ a T ^g2-2(C + Dq 2 )(C + 3Dq 2 )q 2 

d ~ q 2 (C + Dq 2 ) ■ (15) 

The relation (|13|) plays a role of the boundary condition for the equation i|l(J|) for the amplitude. It has a form 
similar to the mentioned in the Introduction boundary condition of de Gennes. Real-valuedness of <5 implies that 
a 2 — 2 (C + Dq 2 ) (C + 3Dq 2 ) q 2 > 0. Taking for the q value the expression JSJ), we obtain 



5 =- A S{ a± r 2 -%)- (16) 

The discriminant a 2 — is always positive since C < 0, D > 0; therefore, S is always real. The two roots of l|14|l 

have opposite signes, but only one of them has a physical meaning at a given magnitude of qL. Indeed, compatibility 
of equations (|11(1 and 112|l implies that 

signS = signt&n-t^-. (17) 
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This consideration suggests that the magnitude of q is rigidly determined by the condition JSJ. We have in this 
case an oscillating dependence of the amplitude logarithmic derivative at the boundary on the the "sample phase 
thickness" qL and, therefore, on the temperature. A behaviour of the system is ruled really by more complicated 
equations than our ones obtained within the slowly-varying amplitude approach. This means that the real system 
can violate some of our assumtions in the boundary layer in search of the energy minimum. One way is to tune q so 
that the energy is minimum. For thick samples qL is rather large and only a very small correction of q would suffice. 
Another possibility is to change the wave number q only within some boundary layer. Only numerical analysis of 
the fourth-order equation can answer the question: which way is chosen by the system. It is easy to see that under 
the assumptions made the problem for the incommensurate phase transition is reduced to a boundary value problem 
which is very similar to that for the ferroelectric thin film. The surface extrapolation length S in l|13l) describes the 
difference between the surface and the bulk and can be positive or negative. It is natural to suppose that for free thin 
films 5 should be positive. Size effect in thin films is noticeable when the extrapolation length value is comparable 
with the film thickness. Thus, if L is assumed to be much longer than the bulk modulation period 2ir /q, the same is 
valid for 5: S ^> 2n/q. We consider in the next section the analytic solution for the second-order equation. 



V. ANALYTIC SOLUTION 



It is fruitful to analyze the phase portrait of the equation (4), which is presented in Fig. 1. 

One can see two regions of solutions separated by separatrices which are shown by dashed lines in Fig. 1. The 
enclosed curves, which correspond to periodic solutions of (4), lie within the separatrices. The boundary condition 
eq. (7) look on the phase diagram as a straight line with negative slope for positive extrapolation lengths and vice 
versa. It is seen from Fig.l that for positive extrapolation lengths all possible solutions lie within the region enclosed 
by the separatrices. 

The analytic solution of the equation (4) reads 
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where sn(z,k) is the Jacoby elliptic sine function [22], uq — u(x — 0) is the value of u in the film middle, k is the 
elliptic module: 
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and uq is related with S by the following relationship: 
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(tn and dn are the Jacoby elliptic functions [22]). 



The obtained solutions allow us to find a distribution of the order parameter amplitude over the film for some par- 
ticular sets of phenomenological parameters, film thickness, and temperature. They also allow us to find a dependence 
of the phase transition temperature Ti on the film thickness for a particular set of the phenomenological parameters, 
Ti stands for the temperature where uq (and therefore u(x)) vanishes. 

The order parameter amplitude u versus x is depicted in Figs. 2 and 3 for some particular positive as well as 
negative values of 5 for several magnitudes of the temperature. The phenomenological parameters used are listed 
in the figure caption. The general behaviour of u(x) for positive S is similar to the behaviour of the spontaneous 
polarization in ferroelectric thin films, that is u decreases near the film surfaces and its value in the film middle is 
smaller than for the relevant bulk (Fig. 2). 

For negative extrapolation lengths there can be two types of solutions for u{x) in agreement with the phase diagram 
(Fig. 1). The behaviour of the solution shown in Fig. 3 is similar to that discussed for ferroelectric thin films: u 
increases near the film surfaces and its value in the film middle is larger than for the relevant bulk. 
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The amplitude u within the film for the second solution is smaller than in the bulk, while its absolute value increases 
near the surfaces, and the amplitude changes its sign at some x that is equivalent to a change of the incommensurate 
modulation phase. However, calculation of the equilibrium thermodynamic potential of the film for two solutions at 
the negative 5 shows that the first kind of solutions (with u is higher than in the bulk) is energetically more favourable. 

The size dependence of the incommensurate phase transition temperature for positive and negative extrapolation 
lengths is shown in Fig. 4. One can see that Ti decreases with decreasing film thickness at positive 5 and vanishes at 
some critical thickness. Similar dependence was obtained for ferroelectric and antiferroelectric thin films, when S was 
implied to be positive [9,17]. The size dependence of Tj for negative S is opposite. 

VI. CONCLUSIONS 

In conclusion, the phenomenological theory of the incommensurate phase transition in thin films is presented. It 
was shown that under the assumption of slowly-varying amplitude of the order parameter the model can be reduced 
to a boundary value problem for the second-order differential equation with conventional homogeneous boundary 
conditions that can be deduced from ones for the fourth-order equation. The amplitude of the order parameter was 
found to decrease near the film surface for positive extrapolation lengths and to increase for negative extrapolation 
lengths. The temperature T"j of the transition was shown to decrease with decreasing the film thickness for positive 
extrapolation lengths, the phase transition cannot occur in films with thickness less than some critical value. 
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FIG. 1: The phase portrait for eq. 4. A=1000, B=100, C=200, D-l, T =l. The dashed lines are the separatrices. The straight 
lines correspond to the boundary condition (7) at 5=1 (1) and 5=-l (2). 
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FIG. 2: The order parameter amplitude u versus x (0 < x < L/2) for positive 5 at various temperatures. L = 2; other parameters 
are given in the caption to Fig. 1. The dashed straight lines show the amplitude in the bulk at the same temperatures. 
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FIG. 3: The order parameter amplitude u versus x (0 < x < L/2) for negative S at various temperatures. L — 2; other 
parameters are given in the caption to Fig. 1. The dashed straight lines show the amplitude in the bulk at the same 
temperatures. 




FIG. 4: Size dependence of the transition temperature. 



